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ABSTRACT 

We present a comprehensive study of X-ray emission and wind properties of massive 
magnetic early B-type stars. Dedicated XMM-Newton observations were obtained for three 
early type B-type stars, CMa, V2052Oph, and £Cas with recently discovered magnetic 
fields. We report the first detection of X-ray emission from V2052Oph and CCas. The latter 
is one the softest X-ray sources among early type stars, while the former is one of the X- 
ray faintest. The observations show that the X-ray spectra of our program stars are quite soft 
with the bulk of X-ray emitting material having a temperature of about 1 MK. We compile the 
complete sample of early B-type stars with detected magnetic fields to date and existing X-ray 
measurements, in order to study whether the X-ray emission can be used as a general proxy 
for stellar magnetism. We find that the X-ray properties of early massive B-type magnetic stars 
are diverse, and that hard and strong X-ray emission does not necessarily correlate with the 
presence of a magnetic field corroborating similar conclusions reached earlier for the classical 
chemically peculiar magnetic Bp-Ap stars. 

We analyze the UV spectra of five non-supergiant B stars with magnetic fields (r Sco, 
/3Cep, i^ 1 CMa, V2052Oph, and £Cas) by means of non-LTE iron-blanketed model atmo- 
spheres. The latter are calculated with the Potsdam Wolf-Rayet (PoWR) code, which treats 
the photosphere as well as the the wind, and also accounts for X-rays. With the exception of 
r Sco, this is the first analysis of these stars by means of stellar wind models. Our models ac- 
curately fit the stellar photospheric spectra in the optical and the U V. The parameters of X-ray 
emission, temperature and flux are included in the model in accordance with observations. We 
confirm the earlier findings that the filling factors of X-ray emitting material are very high. 

Our analysis reveals that the magnetic early type B stars studied here have weak winds 
with velocities not significantly exceeding v csc . The mass-loss rates inferred from the analysis 
of U V lines are significantly lower than predicted by hydrody namically consistent models. We 
find that, although the X-rays strongly affect the ionisation structure of the wind, this effect 
is not sufficient in reducing the total radiative acceleration. When the X-rays are accounted 
for at the intensity and temperatures observed, there is still sufficient radiative acceleration to 
drive stronger mass-loss than we empirically infer from the UV spectral lines. 

Key words: stars: massive - stars: magnetic field - stars: mass-loss - X-rays: stars - tech- 
niques: spectroscopic - stars: individual: t^ 1 CMa, r Sco, (3 Cep, V2052 Oph, £ Cas 



1 INTRODUCTION 

Magnetic fields in massive stars (M*>8Mq) have the potential t o 
influence stellar formation and evolution (e.g. lFerrario et al .120091) . 
Although discoveries of massive star magnetism are increasing 
(Hubrig et al. 2011; Petit et al. 2011; Wade et al. 2011), our knowl- 
edge of their existence, origin and the role they play remains lim- 
ited. The routine presence of magnetic fields in these stars is obser- 
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vationally neither established nor disproved, due to the lack hitherto 
of suitable diagnostic t ools and observations of ad equate quality 
(for a recent review see iDonati & Lan dstreet 2009)). The Zeeman 
effect, commonly used to measure magnetic fields in solar-type 
stars, is less useful in hot stars: the line broadening by macroscopic 
motions exceeds the Zeeman splitting, unless the magnetic field is 
extremely strong. 

While direct measurements of magnetic field have been pos- 
sible only for the closest and brightest stars, there has been indi- 
rect evidence for magnetic fields on massive stars. The wide range 
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of observational phenomena, such as UV wind-line periodic vari- 
ability, cyclic variability in Ha and He II A4686, and excess emis- 
sion in UV-wind lines centered about the rest wavelength (e.g. 
iHenrichs et al. 2005; Schnerr et al .1120081) are commonly explained 
by the influence that magnetic fields exert on stellar winds. For ex- 
ample, cyclical wind line variability is likely due to wind flow that 
is guided by a large scale, dipole-like magnetic field corotating with 
the star. In these cases the timescale of the variability is similar to 
the rotation period. Other indicators of magnetic fields may include 
chemical peculiarity, certain pulsation behavior, non-thermal radio 
emission, and anomalous X-ray emission. 

In this paper we address the X-ray emission and wind prop- 
erties of magnetic B-type stars. Our study includes stars earlier 
than B2 that have both confirmed magnetic fields and X-ray data. 
Thus, our sample includes only massive stars, with initial masses 
exceeding ~ 8 Mq. We consider the B-type stars, where magnetic 
fields have been detected only recently as well as the earliest types 
of "classical" chemically peculiar Bp type stars. The latter hav e 
high incidence of strong magnetic fields (e.g. iBvchkov et al.l2003l) . 
Much work has already been done to study the X-ray emission from 
chemically peculiar magnetic Bp and Ap stars. Using the results 
from X-ray surveys, correlations between the X-ray properties of 
these stars and their ma gnetism were investigated by Drak e et ail 
( 1 19941) and lLeonei l ll994l) . Those works considered a broader range 
of spectral classes, not only B0-B2 stars as this work. The con- 
clusion was reached that, despite some notable exceptions, there is 
little evidence that the X-ray emission from the Bp stars is differ- 
ent fro m other B-stars with similar spectral types. Iczesla & Schmittl 
d2007t) carried out Chandra high-angular resolution observations of 
a sample of late B-type and A-type stars with measured magnetic 
fields in the range from 0.2-17 kG. They showed that the existence 
of a magnetic field of kG strength on a late B-type or A-type star 
is not necessarily a prerequisite for finding X-ray emission among 
these stars. 

The main focus of this work is on j3 Cephei-type variables and 
a slowly pulsating B-type star, where magnetic fields have been 
detected only recently. We performed dedicated observations with 
XMM-Newton for three magnetic B-stars: £ x CMa, V2052 Oph, and 
£ Cas. Two of them, V2052 Oph and £ Cas are detected in X-rays 
for the first time. We also searched the X-ray archives to collect X- 
ray data for other magnetic early B-type stars. Thus, at the time of 
writing, we have assembled a complete sample of early-type mag- 
netic B-stars for which X-ray observations are available. This sam- 
ple will, undoubtedly, grow in future. 

To infer stellar and wind parameters for the five stars in our 
sam ple we employ the non-LTE stellar atmosphere model PoWR 
(e.g. lHamann & Grafeneril2003h . 

Our sample of B-stars stars with directly measured magnetic 
fields is biased to non-supergiant stars. Our program stars allow us 
to explore the parameter space comprising different wind momenta, 
field strengths, and rotation periods with regard to the observed 
X-ray emission. We compare the properties of X-ray emission be- 
tween our program stars and other B-type stars of similar spectral 
types, but without confirmed magnetic field. This allows investiga- 
tion of whether the X-ray emission can serve as an indicator of the 
surface stellar magnetic field. 

In the absence of a magnetic field, the mechanism that pro- 
duces X-ray emission in early-type stars is thought to be the line- 
driven instability (LDI) that is an intrinsic property of stellar winds 
dLucv & Solomoril970h . The LDI mechanism generates numerous 
shocks in a stellar wind, w here plasma can be hea ted up to X-ray 
emitting temperatures (e.g. lFeldmeier et alii 1 997al lbh . 



In addition to the LDI mechanism, X-rays are thought to arise 
when magnetic fields are capable of governing the wind streams. 
Knowledge of the wind parameters, magnetic field topology and 
strength, together with X-ray observations provides a good basis to 
review critically the models of X-ray generation in massive mag- 
netic stars. 

iGroote & Hunger! dl982l) proposed a scenario to explain ob- 
servational phenomena associated with the B2Vp aOriE. Their 
model predicts the formation of a torus of matter around the star, 
formed by a weak stellar wind channeled along magnetic field lines. 
The torus is magnetically coupled to the star. Its co-rotation ex- 
plains the observed photometric and spectroscopic variability. The 
X-ray emission is produced when the wind streams from oppos- 
ing magnetic hemi spheres collide l eading to strong shock heat- 
ing of the plasma. iTownsend et al .1 d2007h further developed the 
IGroote & HungeJ dl982l) scenario by using semi-analytical models 
of the rigidly rotating magnetosphere (RRM) in the limit of very 
strong (~few kG) magnetic fields. Their model predict very hard 
X-ra y emission with temperatures up to 100 MK. 

Babel & Montmerlel dl997bl) envisage a star with a dipole 



magnetic field that is sufficiently strong to confine the stellar wind. 
Collision between the wind components from the two hemispheres 
in the closed magnetosphere leads to a strong shock. The X-ray 
emission from the Ap-type star IQ Aur was explained in the frame- 
work of this "magnetically confined wind shock model" (MCWS), 
and the presenc e of a magnetic field in th e O-type star 8 1 OriC 
was postulated dBabel & Montmerlelll997aT). The direct co nfirma- 
tion of the magnetic field in this star by Donati et al.l d2002l) proved 
that X-rays have diagnostic potential in selecting massive stars with 
surface magnetic fields. Numerical MHD sim ulations in the frame- 
work o f the MCWS model were performed bv lud-Doula & Owockil 
d2002l) . These si mulations compare well with the X-ray observa- 
tions of 8 1 OriC dGagne et al .112005b : however, the model has dif- 
ficulties in describing the X -ray emission from other magnetic O- 
type stars dNaze et al.ll2010l) . In this paper we discuss the applica- 
bility of the MCWS model to the early-type magnetic B stars. 

ICassinelli et alj d2002l). iMaheswararJ d2003l) . iBrown et al] 
d2008l) . and iMaheswaran & Cassinellil d200si) studied the case of 
fast rotating magnetic massive stars with an increasing degree of 
model sophistication to address the formation of disks in classical 
Be-type stars. They showed that magnetic torquing and channel- 
ing of wind flow from intermediate latitudes on a B star could, for 
plausible field strengths, create a dense disk a few stellar radii in 
extent in which the velocity is azimuthal and of order the local Ke- 
plerian speed. Unlike all the other magnetical ly-guided wind sce- 
narios mentioned, the disk model proposed bv lBrown et al.l d2008t) 
is in a steady state, wind inflow being offset by a very slow outflow 
across reconnecting field lines. iLi et al] d2008l) proposed a model, 
where the X-rays are produced by wind material that e nters the 
shock s above and below the disk region. The model by ILi et al] 
d2008l) predicts a relation between the X-ray luminosity normal- 
ized to the stellar bolometric luminosity (LxALboi) and the mag- 
netic field strength in Be-type stars. Compared to the Be-type stars, 
the stars in our sample are slower rotators and have lower wind 
density. 

The paper is organized as follows. Our program stars are intro- 
duced in Sect.|2] The new XMM-Newton observations are described 
in Sect.[3] and analyses of X-ray spectra presented in Sect.|4] The 
X-ray properties of early-type Bp stars are discussed in Sect.[5] 
Analysis of UV spectra for our program stars is presented in Sect.[6] 
We discuss our results in Section|7] with a summary of the main 
points given in Section[8] 
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Table 1. Magnetic early B-type stars with available X-ray observations 



Name 


HD 


Sp 


£ a 
G 


v sin i 
kms 


Prot 

d 


Dipole 


Obliquity 13 




Ref 


T SCO 


149438 


B0V 


(~ 500> 


5 


41.033 


no 




1 


/3 Cep-type and SPB-type stars 


^CMa 


46328 


B0.7IV 


5300 ± 1100 


9±2 


2.18 


yes 


79.1 


2 


/3Cep 


205021 


B2III 


360 ± 40 


27 ±2 


12.00089 


yes 


85° ± 10° 


3, 4 


V2052Oph 


163472 


B1V 


250 ± 190 


60 ±4 


3.63883 


yes 


35° ± 18° 


5 


CCas 


3360 


B2IV 


,,[-+120 


17 ±3 


5.37045 


yes 


77° ± 6° 


6 


Peculiar B stars 


NUOri 


37061 


B0V(n) 


~620 


225 ± 50 




yes 




7 


V1046Ori 


37017 


B2V 


(~ 1500} 


<95 


0.9 


7 


42°-59° 


8, 9, 13 


HR 3089 


64740 


B1.5Vp 


(572 ± 114} 


160 


1.33 


9 




9, 13, 14 


LP Ori 


36982 


B2Vp 


~1100 


80 ± 20 




yes 




7 


a OriE 


37479 


B2Vp 


~ 10000 


140 ± 10 


1.191 


yes 


66° 


10 


HR 5907 


142184 


B2.5Ve 


~20000 


280 


0.5083 


yes(?) 


4° (?) 


11, 12 



a for dipole field configuration, B is the polar field strength. 

For t Sco, HR 3089, and V1046 Ori, an approximate average field strength is shown in angle brackets 
b j3 is the angle between the magnetic and rotational axes for the dipole field configuration 

References: 1 Donati et al. (2006) ; Sota et al. (2011) ; 2 Hubrig et al. (2011) ; 3 Telting et al. (1997) ; 4 Donati et al. (2001) ; 

5 Neiner et al. (2003b) 6 Neiner et al. (2003a) ; 7 Petit et al. (2008) ; 8 Bohlender et al. (1987) ; 

9 Romanvuk & Kudrvavtsev (2008) ; 10 Reiners et al. (2000) ; 11 Abtetal. (2002) ; 12 Grunhut et al. (2010) ; 

13 Bvchkov et al. (2003) ; 14 Borra & Landsfreet (1979) 



2 THE PROGRAM STARS 

Our sample of early type magnetic B-stars with recently discovered 
magnetic fields includes the B0 type star r Sco, the f3 Cephei type 
variables /3Cep, £ CMa, and V2052Oph, and the slow-pulsating 
B-type star (SPB) £ Cas (see TableQ). 

rSco is a well studied object, which has been observed in 
X-rays by all major X-ray missions. A complex magnet ic field 
topology was discovered in rSco by IDonati et all feOOfjh . r Sco 
is a source of hard X-ray emission and displays n arrow emis- 
sion l ine profiles, suggesting nearly stationary plasma dMewe et al .1 
120031) . Recently, rSco was monitored throughout its rotational 
cycle by the Suzaku X-ray observatory. Contrary to expecta- 
tions, rotational mo dulations of X-ray emission were not detected 
dlgnaceetai]|2010h . Magnetic fields were detected on two other 
stars, H D 66665 and HD 63425, that are spectroscopically similar 
to t Sco dPetit et al ■1201 ll) . Although it remains to compare their X- 
ray characteristics to r Sco, the new discoveries suggest that r Sco 
may be a prototype for a wider class of stars. 

/3Cep, £ x CMa and V2052Oph are [3 Cephei type variables. 
ft Cep-type stars have masses between 7 and 16 Mq and are on 
the main-sequence or in an early post-MS evolutionary phase 
(corresponding to spectral type 09 to B3). Stars of this type 
display photometric and spectral variations caused by pulsations 
in low-order pressure and gravity modes of short-period (3-8 h) 
dDziembowski & Pamiatnvkhll 19931) . 

A magnetic field in f3 Cep was detected by iHenrichs et al.l 

d2000h. Based on the analysis of the Rosat measurements, 
IDonati et all ( 1200 ll) suggested that the X-rays from ft Cep can be 
explained within the framework of the MCWS model. This model 
predicts rotationally modulated, strong and hard X-ray emission 
originating at a few stellar radii from the surface close to the disk 
around the magnetic equator. However, data appear not to support 
these predictions. X-ray multiphase high-spectral resolution obser- 
vations of f3 Cep by XMM-Newton and Chandra fail to show ei- 
ther rotational modulations or modulations on the time-scale of the 



stellar pulsations dFavata et al.ll2009h . The X-rays do not originate 
further out in the wind than 2-3 7?*, and the emission is quite soft 
with the bulk of X-ray emitting plasma at T«3 MK, dFavata et alj 
l2009h . The X-ray luminosity of (3 Cep is not atypical of stars with 
the same spectral type. 

iHubrig et all d2006|) and lSilvester et all d2009l) measured mag- 
netic fields in a sample of eight f3 Cep-type stars. Both these studies 
agree in the lack of a definite detection of magnetic fields in the six 
f3 Cep-type stars, with typical longitudal field formal errors of few 
tens of G. The results of these investigations seems to indicate that 
the presence of a strong magnetic field is not a general intrinsic 
property of this spectral class. 

Magnetic field measurements of CMa were report ed in 
IHubrig et all d2006l) and c onfirmed by ISilvester et ail d2009h . Re- 
cently, I Hubrig et al .11201 ll) established the magnetic field configu- 
ration, and the rotational periods for a number of magnetic B-stars, 
among them £ CMa. It was shown that we see £ CMa nearly ro- 
tational pole-on. 

j 1 CMa was prev i ously observed in X-rays by Einstein and 
Rosat. ICassinelli etail d 19941) presented a Rosat PSPC spectrum 
of this star. They reported a temperature of 3.7 MK and an emis- 
sion measure EM~ 10 54 cm~ 3 , and pointed out that £ x CMa has 
the highest X-ray luminosity and the hardest spectrum among the 
[3 Cep-type stars in their sample. 

The detection of a magnetic field in V2052 Oph was reported 
bv lNeiner etail d2003bl) . The polar magnetic field was determined 
assuming that the star is an oblique rotator. iNeiner et"aT 1 (l2003bl) 
note that V2052 Oph is very similar to (3 Cep. They suggested that 
based on f3 Cep, the wind of V2052 Oph should be confined and 
form a disk-like structure. Prior to our observations, there was no 
positive detection of X-ray emission from V2052 Oph, with the 
Rosat AH Sky Surx'ey (RASS) yielding only an upper limit to its 
X-ray luminosity. 

The slow-pulsating B stars are slightly less luminous and 
cooler than the {3 Cep-type stars. These stars show multi-periodic 
brightness and color variations on a time-scale of 0.8d-3d. 
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Table 2. XMM-Newton observations of three magnetic early B-type stars 



Star 


MJD 


useful exposure 


PN count-rate a 






[ksec] 


[s- 1 ] 


^CMa 


55046.4576 


6.9 


0.61 ± 0.01 


V2052Oph 


55077.0906 


9.2 


0.003 ± 0.001 


£Cas 


55046.4576 


12.3 


0.040 ± 0.002 



in the 0.3-7.0 keV band; background subtracted 



iNeiner et al.l d2003al) reported the disco very of a magnetic field on 
the SPB star C Cas. iHubrig et alj feOOfj) and lSilvester et al.l d2009t) 
searched for magnetic fields in a large sample of SPB stars; how- 
ever, the results of their measurements agree only partially - both 
studie s detect a magnetic field on 16 Peg (B3V) ^Silvester et al.l 
d2009h report a marginal detection). 16 Peg is not yet detected in 
X-rays. Prior to our observations, there were no positive detection 
of X-ray emission from £ Cas eit her, with the RASS y ielding an up- 
per limit to its X-ray luminosity. INeiner et al.l d2003al) analyzed the 
magnetic field configuration of £Cas and concluded that the star 
is an oblique magnetic dipole. However, they argued that there is 
no evidence for a disk around this star, based on the lack of an IR 
excess. 



3 OBSERVATIONS 

We obtained dedicated XMM-Newton observations of £ CMa, 
(Cas, and V2052Oph. All three (MOS1, MOS2, and PN) Euro- 
pean Photon Imaging Cameras (EPICs) were o perated in the stan- 
dard, full-frame m ode and a thick UV filter dTurner et al.l 1200 ll ; 
IStriider et al. I l200lh . The log of observations is shown in Table|2] 
The data were analyzed using the software SAS 9.0.0. Each of the 
stars in our sample was detected by the standard source detection 
software. The exposure times and EPIC PN count rates for our pro- 
gram stars are given in Table[2] The spectra and light-curves were 
extracted using standard procedures. The spectra of all sources 
were extracted from regions with diameter ~ 15". The background 
areas were chosen in nearby areas free of X-ray sources. The EPIC 
PN spectra of £ x CMa were corrected for "out-of-time events". 



4 X-RAY SPECTRA AND THEIR FITTING 

To anal yze the spectra we used the standard spectral fitting software 
XSPEC dArna ud 1996) . The referen c e abun dances were set to solar 
values according to lAsplund et al. I d2009l) . The number of counts 
per bin in the spectra o f C Cas and V2052 Oph is small, therefore 
we used Cash-statistic dCashlll979l) to fit their spectra, while the 
£ J CMa spectra were fitted using the x 2 -statistics. Using the neu- 
tral hydrogen column density as a fitting parameter does not yield 
a sensible constraint on its value. Therefore, Nn was fixed at its in- 
terstellar value as found from the analysis of optical and UV spectra 
(see Section[6]l for all our program stars. 

The standard model vapec, which assumes that the plasma 
is in collisional ionization equilibrium (CIE), was used to model 
the observed X-ray spectra of our program stars. CIE models 
require that several equilibria exist in the stellar wind namely 
that the electrons and ions are both thermalized and have equi- 
librated so as to have T c = 71 and that ionization and recom- 
bination rates are balanced. The thermalization processes is con- 



trolled by Coulomb collisions or even faster plasma effects. How- 
ever, cooling of the plasma can be faster than the overall recom- 
bination time due to of som e ions with important cooling lines. 
ISutherland & Dopital d 19931) provide a comparison of the colli- 
sional equilibrium timescale and the cooling timescale for carbon 
and iron ions. From their Fig. 15, at temperatures above > 1 MK, the 
cooling time is r coo i = l.5nkT / 'n B n\h(T) > 5 x 10 8 /n e ) hr, and 
the collisional equilibrium time for C VI is r rcc0 mb = a re J omb • 
n^ 1 < 1.5 x 10 8 (1 cm~ 3 /n c ) hr, where n c and ni are elec- 
tron and ion number densities respectively in cm~ 3 , A(T) is the 
cooling function, and a rGCOm b is the recombination coefficient (see 
ISutherland & Do pita ( 1993) for details). Hence, at the temperatures 
where the maximum of emission measure is determined from the 
analysis of X-ray spectra (see Table|3}, the cooling of the plasma 
is slower than the recombination time-scale, and the conditions for 
the establishing the CIE seem to be fulfilled. 

Using CMa as an example, we deduce from our wind mod- 
els (see Section[6} that at 1 R* distance from the photosphere, the 
electron density in the wind is ~ 10 7 cm -3 . The cooling time of a 
1 MK plasma is ~ 50 hr, while the collisional equilibrium time for 
C VI is ~ 15hr. Note that this is long compared to the pulsational 
time-scale of /3C ep-type variables, e.g. the pulsational period of 
£ 1 CMa is 5.03 hr dStankov & Handlerll2005l) . 

The mean free path for a Coulomb collision Z m f p ~ 9.4 x 



10 7 T 2 



cm, where T is in the units of MK, and n$ is in the 



units of 10 s cm -3 . Using parameters for £ CMa close to the stellar 
surface, the value of i m f p is only 100 cm. Therefore, we believe that 
the use of CIE models to describe the X-ray spectra of our program 
stars is justified. 

The chemical abundances of our program stars deviate from 
the solar, which is most likely a consequence of the surface mag- 
netic field. Peculiar abundances are often found in magnetic stars, 
and is commonly explained by diffusion processes which allow 
heavier elements to sink in the atmosphere under the influence of 
gravity, while lighter elements are lifted to the surface by radiation 
pressure. The summary of the relevant publications with analyses 
of abundances in our program st ars can be found in INeiner et al.l 
d2003bl) , lNeiner et alJ l EoOGah and lMorel et"al1d2008l) . In summary, 
these studies agree that the magnetic stars show an overabundance 
of nitrogen, and sometimes helium. The quality of the X-ray spectra 
are not sufficient to constrain abundances. T herefore, for each star 
we have used the abundances determined bv lMorel et all d2008l) . 

^ CMa The observed EPIC spectra of CMa and the fitted model 
are shown in Fig.Q] The abundances were set to solar v alues, except 
C/C a = 0.56, N/N =1.48, and O/O =O.79, as found bv lMorel et ail 
d2008l) . We initially at tempted to fit the sp ectra using two tempera- 
ture models similar to iFavata et alj d2009l) for modeling the /3 Cep 
X-ray spectrum. However, we found that adding a softer model 
component with kT w 0.1 keV significantly improved the qual- 
ity of the fit. Parameters of the best-fit model are shown in Ta- 
ble^ We also attempted to fit the £ x CMa spectrum assuming a 
higher temp erature component . By analogy with the r Sco spec- 
tral models dMewe et alj|2003l) . a spectral model component with 
kTi = 1.7keV was assumed. However, the emission measure of 
this hot component cannot be constrained. As a next step, we in- 
vestigated whether the £ x CMa spectrum can be fitted with non- 
equilibrium ionization models. Using the nei model available in 
XSPEC, we failed to find a suitable set of model parameters to re- 
produce the observed spectrum. 

V2052Oph Our XMM-Newton observation detected the X-ray 
emission from V2052 Oph for the first time. The source has a very 
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Energy (keV) 

Figure 1. XMM-Newton PN (upper curve), and M0S1 and MOS2 (lower 
curves) spectra of CMa with the best fit three-temperature model (solid 
lines). The model parameters are shown in Table|3] 
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Figure 2. XMM-Newton EPIC-PN spectrum of V2052 Oph and the best fit 
two-temperature model. The model parameters are shown in TablefJ] 



low count rate (see Table[2}; in total only about 60 source counts 
were collected. At a distance of d = 417 pc, the X-ray luminosity 
of V2052 Oph is only L x ~ 3x 10 29 ergs -1 , making it one of the 
least X-ray luminous early type B-stars in the sky. To model the 
observed low S/N spectrum, we fixed the neutral hydrogen column 
density at its interstellar value. The PN spectrum of V2052 Oph 
can be well described using a two temperature plasma model see 
Tabled. The abunda nces C/C a =0 6, N /N =1.4, and O/O Q =0.5 
were used as found bv lMorel et al.l J2008h . The temperature and the 
emission measure of soft component are only poorly constrained 
but the presence of a soft component is required to reproduce 
the observed hardness ratio. There are no indications of a harder 
component being present in the spectrum of V2052 Oph. 

Pausing to consider the X-ray spectral analyses of the /3 Cep- 
type stars (/3Cep, ^ CMa and V2052Oph), it emerges that the 
spectral temperatures are similar among these stars. We define the 
mean spectral temperature (kT) as the emission measure weighted 
average temperature, with 

(kT) = Y^kT r - EM,/ EM, (1) 

i i 

Then, the mean spectral temperature for all magnetic p Cep-type 
stars is the same at ~3.5 MK (see Table[3]). However, we find a large 
difference in X-ray luminosities for the /3 Cep stars, with £ CMa 
being 4 times more X-ray luminous than P Cep and 50 times more 
X-ray luminous than V2052 Oph (see TableQJ. 

£ Cas Our XMM-Newton observations of £ Cas are the first X-ray 
observations to yield a spectrum of a magnetic SPB-star. Visual in- 
spection of EPIC XMM-Newton images already reveals that £ Cas 
is a very soft X-ray source. This expectation is confirmed by the X- 
ray spectral fits. The EPIC spectra of £ Cas and the best fit 2T CIE 
model are shown in Fig.|5] A solar compositio n was assumed, ex - 
cept for C/C Q =0.5, N/N =1.4, and O/O Q =0.5 jMorel et alj|2008h . 
Parameters of the model are shown in Table|3] The maximum tem- 
perature inferred from the spectral analysis is ~4 MK, and there are 
no indications of a harder spectral component. The emission mea- 
sure is dominated by the 1 MK plasma; the hotter, 4 MK component 
constitutes less than 20% of the total emission measure. £ Cas has 
the softest X-ray spectrum among all hot stars where magnetic field 




Energy (keV) 

Figure 3. XMM-Newton EPIC PN (upper), and MOS1 and MOS2 (lower) 
spectra of £ Cas and the best fit two-temperature model. The model param- 
eters are shown in TablefJ] 

have been detected. The mean X-ray spectral temperature of £ Cas 
is about 1 MK (see Tabled. 



5 X-RAY PROPERTIES OF CLASSICAL PECULIAR 
MAGNETIC EARLY-TYPE B STARS 

In this section we briefly consider the X-ray properties of chem- 
ically peculiar Bp-stars, addressing only stars with early spectral 
types. The classical chemically peculiar Bp-Ap stars have a high 
incidence of strong surface magnetic fields. The observed temporal 
variations of longitudal magnetic fields are usually compatible with 
dipole or low-order multip ole fields inclined to the rotation axis 
(Donati & Landstreetl 12003) . The strong magnetic fields strongly 
influence weak stellar winds, making the Bp stars excellent test 
cases for models of X-ray production in magnetic early-type stars. 

There has been earlier work se eking to establish the X-ray 
properties of chemically peculiar stars. iDrake et al.ldl994h searched 
the RASS database at the positions of about 100 magnetic Bp-Ap 
stars. They detected 10 X-ray sources and argued that in four cases 
the X-ray emission presumably arises from an early-type star with 
a radiatively driven wind. The X-ray luminosities were found to 
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Table 3. The spectral parameters derived from the XMM-Newton EPIC observations of our program stars assuming the multi-temperature CIE plasma models 
(vapec) corrected for the interstellar absorption (tbabs). The values which have no error have been frozen during the fitting process. The corresponding spectral 
fits are shown in Figs. |l|3|2l For comparison, the spectral parameters inferred from modeling XMM-Newton data of /3 Cep and Suzaku spectra of r Sco are also 
shown. 



Star 


£ J CMa 


CCas 


V2052 Oph 


/3Cep a 


rSco b 


[10 20 era" 2 ] 


1.4 


3 


15 


2.5 ±0.1 


3 


hTi [keV] 


0.12 ± 0.01 


0.08 ± 0.02 


0.14 ±0.12 


0.24 ±0.01 


0.11 ±0.01 


EMi [10 53 cm- 6 ] 


22.48 ± 5.64 


1.29 ± 1.05 


0.006 ± 0.025 


11 ± 2 


17.0 ±2.61 


kT 2 [keV] 


0.32 ± 0.01 


0.31 ± 0.02 






0.34 ±0.01 


EM 2 [10 53 cm- 3 ] 


19.3 ± 3.36 


0.27 ±0.07 






10.4 ±0.51 


kT s [keV] 


0.68 ±0.05 




0.65 ± 0.11 


0.69 ±0.03 


0.71 ±0.10 


EM 3 [10 53 cm" 6 ] 


6.41 ± 2.57 




0.003 ± 0.002 


1.3 ±0.3 


7.2 ±0.3 


kTi [keV] 










1.52 ±0.06 


EM 4 [10 53 cm" 6 ] 










5.2 ±0.3 


(kT) ee Ei kT t ■ EM,/ E» EM, [keV] 


0.3 


0.1 


0.3 


0.3 


0.5 


Flux d [10~ 12 ergcm~ 2 s" 1 ] 


1.1 


0.093 


0.006 


1.0 


16.4 



a the values are adopted from Favata et al. (2009) 

b the values are adopted from Ignace et al. (2010) 

c correspond to the ISM hydrogen column density for all stars 

d in the 0.3-7.0 keV band, except of r Sco in the 0.3-10.0 keV band, absorbed; 



be in general agreement with the canonical values of X-ray emis- 
sion from massive stars Lx ~ 10~ 7 Lboi- E.g., the B1.5Vp star 
HR3089 has logL x /L bol w -7.2. 

To include newer data, we selected Bp stars with spectral 
types earlier than B2 from The Ca talog of Peculiar Magnetic Stars 
dRomanvuk & Kudrvavtsevl|2008i ). This sam ple was augmente d by 
two more magnetic Bp-type stars reported bv lPetit et al.l 1( 120081) and 
the magnet ic Be star HR 5907 re cently reported by the MiMeS col- 
laboration dGrunhut etaill2010h . As a next step we searched the 
available X-ray catalogs. Our search revealed that a OriE, LPOri, 
NU Ori, and V1046 Ori, HR 5907 have positive X-ray detections. 

a Ori E is the prototypical oblique fast magnetic rotator. Its X- 
ray emission is quite hard: the EM of 10MK plasma is as large as 
the E M of the cooler plasma with T~ 3 MK dSanz-Forcada et al] 
120041) . The mean spectral temperature of a Ori E (in quiescence) 
is 0.66 keV or 7.6 MK. We estimate the bolometric luminosity of 
<j OriE, log'Lboi/^o = 3.8 , based on its sp ectral type and UBV 
colours and assuming a distance of 640 pc dHunger et al.lll989l) . 
The average (quiescent+flare) X-ray luminosity in 0. l-2.4keV 
band at 640 pc is Lx ~ 8 x 10 31 ergs -1 , the quiescent X-ray lu- 
mino sity is factor of 5 lower than in the peak dSanz-Forcada et all 
120041) . The star is X-ray luminous: its ratio of X-ray to bolomet- 
ric luminosities Lx/iboi is nearly two orders of magnitude higher 
than normally found in B-type stars (see Table|4j- An extremely 
rare event in the history of X-ray observations of massive hot stars - 
an X-ray flare - has been observ ed from a Ori E dGroote & Schmitj 
120041 : ISanz-Forca da et al. 2 0041) that may be underst ood as a result 
of centrifugal breakout of the torus dMullanll2009l) . These X-ray 
propert ies of a Ori E seems to confirm the expectations of the RRM 
model dTownsend et all 20071) . 

Given the general nature of these models, one would naturally 
expect that these results could be applied to other BV stars with 
strong magnetic fields. This however, seems not to be the case. 

The young star LP Ori is similar to o Ori E in many respects: 
it has a similar age, spectral type, and a kG strong magnetic field. It 
is also a source of hard X-ray emission dGetman et al .120 05). How- 



ever, the X-ray luminosity of LP Ori is significantly lower than the 
X-ray luminosity of a Ori E. In fact, the X-ray luminosity of LP Ori 
is nearly the lowest one among all magnetic early B-stars (see Ta- 
ble|4j. The models predicts modulations of the X-ray luminosity 
due to the occultation of X-ray emitting site by the opaque torus. 
These modulations are thought to occur on a time-scale comparable 
to the rotati onal period of the sta r. LP Ori was observed by Chandra 
for «8.8 d dGetman et ai1l2005h . This long exposure time suggests 
that the observed X-ray flux is somewhat averaged over the rota- 
tional period, and the star is indeed intrinsically faint in X-rays. 

X-ray observations of a OriE and LPOri seemingly suggest 
that the strong magnetism may account for the hardness of the X- 
ray emission. However, the Bp stars V1046 Ori and NU Ori provide 
counter examples. 

Similar to a Ori E, V1046 Ori has a kG magnetic field. It also 
is an oblique magnetic rotator, whose obliquity and rotational pe- 
riod are similar to that of o Ori E. The star was serendipitously ob- 
served by XMM-Newton. No X-ray spectra are available, but hard- 
ness ratios are provided in The Second XMM-Newto n Serendipitous 
Source Catalog or 2XMMi-DR3 and in lNazel d2009l) . We have com- 
pared these hardness ratios with those of r Sco reported in the same 
catalog. V1046 Ori is a softer X-ray source compared to r Sco, and 
its X-ray luminosity is relatively low (see Table|4j. We estimate its 
bolometric luminosity based on UBV colors and spectral type. 

The magnetic field on NUOri is only about ~0.5 kG and is 
weake r than the fields of the other Bp stars in our sample. IPetit et all 
d2008l and references therein) comment on NU Ori being a triple 
system, containing a massive B0.5V primary, along with a low 
mass spectroscopic comp anion and a visual comp anion. NUOri 
was recently analyzed by I Simon-Diaz et al.l d201 ll) . who obtained 
logLboi/i© = 4.4. NUOri was observ ed by Chandra simulta- 
neously with LP Ori dGetman et aT1l2005l) . Surprisingly, that anal- 
ysis indicates that the X-ray luminosity of NUOri is a factor of 
30 higher than for LP Ori, yet the former has a considerably softer 
emission as compared to the latter. 

HR 5907 has very strong magnetic field dGrunhut et al .120101) . 
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Table 4. X-ray luminosities of magnetic early B-type stars 



Name d L x log (Lx/Xbol) 

pc 10 30 ergs _1 
tSco 150 40 -6.4 

Magnetic Cep-type and SPB-type stars 
IFcMa 420 30 

/3Cep 200 6.4 -7.0 

V2052Oph 400 0.3 -8.0 

CCas 180 05 -7.5 

Other magnetic early-type B stars 
NUOri 400 ii) ^8 

V1046Ori a 380 0.1 -8.0 

crOriE b 640 80 -5.6 

HR3089 C 300 2 -7.2 

LPOri 470 0.02 -8.5 

HR5907 120 0.4 -7.4 



Distances are from^an Leeuwen (2007) except of LPOri, cr OriE, HR3089; 

a the distance between the X-ray source; 

2XMMJ053522.0-042938 and the position of V1046Ori is 3"; 

b Lx (quiescent+flare) in 0. 1-2.4 keV band from Sanz-Forcada et al. (2004) ; 

c L x in 0.1-2.4 keV band from Drake et al. (1994) ; 

Assuming a dipole, the field strength at the pole is 20 kG. Th e 
star is a fast rotator with vsini « 280 kms -1 jAbt et aljr2002l) . 
lYuditj J200lh detect intrinsic polarization on the level 0.17% in V 
band and consider this star as a classical Be star with a disk. The 
star was observed by Rosat PSPC for « 4.8 ks. It is detected with 
a flux « 2 x 10 -13 erg cm -2 s _1 corresponding to rather low X- 
ray luminosity (see TableHJ. We fit the low signal-to-noise Rosat 
PSPC spectra and find that HR 5907 has a relatively hard spectrum 
characterized by a temperature above 1 keV. 

The rigidly rotating magnetosphere scenario makes robust 
predictions of strong, variable, and hard X-ray emission. It appears 
that in five out of the six early-type magnetized B stars under dis- 
cussion, these predictions are not fulfilled. Considering just our 
small sample of early-type B stars with kG strong magnetic fields, 
we must conclude that the new measurements confirm earlier re- 
sults on the X-ray emission from Bp-Ap stars. The properties of 
their X-ray emission are quite diverse, with majority of the stars 
emitting X-rays on the lev el typical for the general sample of B- 
stars (e.g lDrake et al.lll994l) . 

Placed in the context of this earlier work, our results corro bo- 
rate the basic conclusion that while strong and hard X-ray emission 
is sufficient to suggest a star may be magnetic, it is not a required 
property of magnetic stars. 



6 STELLAR WINDS 

The strong UV radiation of B-type stars drives their stellar winds. 
The parameters of stellar winds are commonly inferred from the 
analyses of optical and UV spectral lines by means of stellar atmo- 
sphere models. While the B-type giants and supergiants have often 
been analyzed, the w inds of BIV a nd BV type stars yet remained 
relatively unexplored. |PrinialJl989h used profile fits to model spec- 
tral lines of Civ, Si IV, and Si III in high-resolution spectra of 
B stars obtained with the International Ultraviolet Explorer (IUE), 
and produced a homogeneous set of wind-velocity and column- 
density measurements for 40 non-emission, non-supergiant B stars. 
It was concluded that the presence of C IV resonance doublet in 
their UV spectra must be due to the effect of X-rays on wind ion- 



Table 5. UV observations of program stars used in the analysis 



object 


data set 


date 


CCas 


SWP51309 


04/07/1994 


C 1 CMa 


SWP 03574 


12/12/1978 


r Sco 


SWP 55997 


23/09/1995 


V2052 Oph 


SWP 50431 


31/03/1994 


/3Cep 


SWP40477 


28/12/1990 



ization . This effect is often referred to as "superionization". IPrinial 
dl989h did not find any asymmetry in the Si IV lines which would 
be typical if the stellar wind was strong. Estimates of column den- 
sities revealed that the mass-loss rates in non-supergiant B stars 
are significantly lower than in O and Be stars. From fitting the UV 
resonance lines of C IV and Si IV, it was found that the products of 
mass-loss rate and ionization fraction, Mq, are <C10 -10 Mq yr _1 . 
Regarding the wind velocities, it was shown that the winds of 
"normal" non-supergiant stars are slow compared to those of the 
Be stars. The maximum observed wind velocities do not generally 
exceed the stellar escape velocity. 

In this paper we analyze the UV spectra of our sample mag- 
netic stars by means of a stellar atmosphere code. Stellar spec- 
troscopy codes used to analyze massive stars are based on spher- 
ically symmetric or plane-parallel geometries, despite of the clear 
observational evidence t hat massive star winds can depart from 
spherical symmetry (e.g. lKaperetal.lll997l ; lHamann et al]|200ll) . 
It is very probable that magnetic fields affect the wind geometry. 
However, at present 3D modeling of stellar winds that would allow 
quantitative spectroscopic analysis is beyond reach. Our approach 
allows to conduct a multiwavelength X-ray and UV/optical spectral 
analysis. The obtained wind parameters of our sample stars can be 
compared to other massive stars in the general framework of mas- 
sive star wind quantitative spectroscopy. 



6.1 Model description: PoWR stellar atmosphere 

In this work we analyze the spectra of magnetic B stars by means 
of stellar atmosphere models that account for the presence of X- 
rays. We use the Potsdam Wolf-Rayet (PoWR) model atmosphere 
code dHamann & Grafenej|2003ll2004 . The PoWR code has been 
used extensively t o analyze not only m assive stars with strong 
stellar winds (e.g. iLiermann et al. but also low-mass cen- 

tral stars of planetary nebulae, and extreme helium stars dHamanrJ 
l20irilTodt et al.l2010l) . The PoWR code solves the non-LTE radia- 
tive transfer in a spherically expanding atmosphere simultaneously 
with the statistical equilibrium equations and accounts at the same 
time for energy conservation. Complex model atoms with hundreds 
of levels and thousands of transitions are taken into account. The 
computations for the present paper include complex atomic mod- 
els for hydrogen, helium, carbon, oxygen, nitrogen, and silicon. 
Iron and iron-group elements with millions of lines are included in 
the P oWR code through the concept of superlevels dGrafener et al.l 
120021) . The extensive inclusion of the iron group elements is impor- 
tant not only because of their blanketing effect on the atmospheric 
structure, but also because the diagnostic wind lines in the UV (e.g. 
the C IV and Si IV resonance lines) are heavily blended with the 
"iron forest". 

The PoWR models can take account of stellar-wind clump- 
ing i n the standard volume-filling factor 'microclumping' approach 
(e.g. lHamann & Koesterkelll998l) . or even, with an approximate 
correction for wind clumps of arbitrary optical depth ("macro- 
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Figure 4. Optical CTIO spectram of CMa obtained on 1988/1 1/02 (Walborn & Fitzpatrick 199ol) (blue line) vs. a PoWR model spectrum (red line) using 
stellar parameters of § 1 CMa as shown in Table[6] The model spectram is convolved with a 2.5 A Gaussian to match the spectral resolution of the "2D-Frutti" 
spectrograph at CTIO. 



clumping". lOskinova et alj|2007l) . However, since the wind signa- 
tures are quite weak for all stars of our sample, and since there 
is no way to estimate the degree of clumping, we apply here only 
smooth-wind models. Note that microclumping does not directly 
influence the resonance line profile shapes via radiative transfer ef- 
fects, but can indirectly affect the lines by virtue of modifying the 
ionization stratification. 

The X-ray emission and its effects on the ionization structure 
of the wind is in cluded in the PoWR atmosphere models according 
to the recipe of lBaumetaLldl992l) . We assume an optically thin, 
hot plasma component distributed within the "cool" stellar wind. 
The uniform value of the so-called Emission Measure 'filling fac- 
tor' , X flu = EMh t I EM coo \ (really a density squared weighted 
volume ratio), is adjusted such that the emergent X-ray luminosity 
agrees with the observations. The X-ray emissivity is restricted to 
the fast-wind domain, for which we assume a minimum radius of 
1.1 R*. The absorption of the X-ray radiation by the "cool" stellar 
wind is taken into account, as well as its effect on the ionization 
stratification by the Auger process. 

The lower boundary of the model atmosphere is set at a Rosse- 
land depth of 20, meaning that the (nearly) static photosphere is in- 
cluded in the computations. The velocity field consists of two parts. 
In the photospheric part of the atmosphere, a hydrostatic density 
stratification is assumed, while for the wind the usual "/3-law" pre- 
scription is adopted, in our case with an exponent @ = 1. We also 
tested the so-called double-/3 law dHillier & MilleJ 19991) . but with- 
out substantial improvement of the fit quality. 

Each stellar atmosphere model is defined by the effective tem- 
perature, surface gravity, luminosity, mass-loss rate, wind terminal 
velocity Vao , and chemical composition. The gravity determines the 
density structure of the stellar atmosphere below and close to the 
sonic point. From the pressure-broadened profiles of photospheric 
lines, the spectroscopic analysis allows derivation of the gravity 
and thus the stellar mass. For our sample of B stars, we confirm a 
discrepancy between the spectroscopi c and evolutionary mass , the 
former being lower than the latter (see lWeidner & VinklfcOlOR and 
references therein). 

Dedicated analyses of the photospheric s pectra o f our program 
stars ha ve b een performed by INeiner et al.l d2003bl) . iNeiner et al .1 
d2003ah and lMorel etail d2008l) . Using their abundances, T eS , and 
log g we find that our model reproduces photospheric spectra very 
well, as demonstrated by the example shown in Fig.|4] Therefore, 
to speed up the analysis, we adopt the literature values of T 8 ff> and 



log p. The effects of rotation with values -usini from Table[T]are 
also included in the models. 

All stars in our sample have measured parallaxes. Hence we 
can scale synthetic spectral energy distributions (SED) to that dis- 
tance and fit to the observations, covering the whole wavelength 
range from the UV to IR. Furthermore, the model spectra are cor- 
rected for interstellar extinction. Dust extinction is tak en into ac- 
count using the reddening law of ICardelli et al.l |l989). An exam- 
ple of a SED fit is shown in Fig.|5] The figure also shows the model 
continuum without lines, illustrating how the "forest" of iron lines 
forms a pseudo-continuum in the UV that is well reproduced by our 
model. 

With logg, T e ff, and Lboi being fixed, we compare the syn- 
thetic and observed lines varying the wind parameters M and Voo 
in order to achieve the best fit. The usual indicators of mass loss 
used in B supergiants and Be stars, as well as in O-type stars, are 
Ha together with the UV resonance lines. However, in our sam- 
ple of non-supergiant stars the Ha line is entirely photospheric^- 
Therefore, our wind diagnostic can only be based on the UV res- 
onance lines. High resolution IUE spectra were retrieved from the 
archive for all stars of our sample. For some of them, multiple IUE 
observations a re available (see e.g . the IUE time series for j3 Cep 
and £ : CMa in lSchnerr et al.ll2008l) . In those cases we selected the 
observation of best quality (cf. Table|5}- 

6.2 rSco 

r Sco is a well studied object that has remained one of the primary 
targe ts in stellar UV and X-ray astronomy from their early d ays 
(e.g. iRogerson & Upsonlll977l : iMacfarlane & Cassinelli|[l989l) . It 
was one of the first stars whose UV spectrum was analyzed by 
means of atmosp here models base d on co-moving radiative trans- 
fer techniques by iHamannl dl98ll) . who obtained a mass-loss rate 
log M — —8.9 ± 0.5 from the analysis of UV lines. It was also 
shown that O VI and N V lines cannot b e reproduced by the cool 
wind models. Icassinelli & Olsonl dl979l) explained the presence of 

VI and N V as arising from Auger ionization by X-rays. 

We re-analyzed the wind properties of r Sco using the P oWR 
model. The abundances as presented in iHubrig et al.l d2008l) were 

1 (3 Cep shows Ha emission episodes. However, ISchnerr et al.l d2006l) 
found that the Ha emission is not related to the primary in Cep, but is 
due to its 3.4 mag fainter companion that is a classical Be star. 
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Figure 5. Spectral energy distribution for £ C Ma. Observed IUE sp ectra are s hown as thin blue l i nes. B lue boxes indicate observed photometry (labels: 
magnitudes) taken from the 2MASS catalogue iSkrutskie et~al]|2006l) and from Morel & Magnenat 1 1978). The synthetic spectrum (red line) is calculated 
using parameters given in Table|6] The model flux is reddened with i?B- V = 0.04 and corrected for interstellar Lyman line absorption. The model continuum 
without lines (red dotted) is also shown to demonstrate how the iron group lines form a pseudo-continuum in the UV. 
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Figure 6. r Sco: The effect of ionization by X-rays on Civ 
AA1548.2, 1550.8 A doublet. Detail of the UV spectrum observed with IUE 
(blue thin line) vs. PoWR models: without X-rays (green dotted line) and 
with X-rays (red thick line). The model parameters: log(M) = —9.3, 
Vao = 1000 km/s. This figure shows that C IV is efficiently destroyed by 
X-rays in the outer parts of the atmosphere. Without accounting for the ion- 
ization by X-rays, the mass-loss rate would be underestimated. 
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Figure 7. r Sco: Same as in Fig.[6]but centered on the N V doublet. 



Figure 8. r Sco: Same as in Figs.|6]and|2]but centered on the Si IV doublet. 
The model parameters: log(Af) = — 8.6, i)oo = 1000 km/s. Note that even 
without X-rays, most silicon is already ionized to Si V in the outer parts of 
the atmosphere, therefore the Si IV doublet is insensitive to X-rays. 
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Figure 9. r Sco: Determination of the wind velocity from modeling the 
C IV line. Detail of the UV spectrum observed with IUE (blue thin line) 
vs. PoWR models with Voo = 500 km s- 1 ( green dotted line) and Voo = 
1000 km s -1 (red dotted line). The derived mass-loss rate is log(M) = 
—9.3. The blue absoiption wing of the C IV line is better matched by the 
model with higher . Both models include ionization due to X-rays with 
parameters as listed in Table [T] 
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Figure 10. r Sco: Constraining the mass-loss rate from modeling of the 
C IV doublet. Detail of the UV spectrum observed with IUE (blue thin line) 
vs. PoWR models with log M = —8.6 (green dotted line) and logM = 
—9.3 (red dotted line). Both models use Voa = 1000 km/s and include 
superionization. The C IV doublet is better matched by the model with the 
lower M. 



adopted. We obtained wind parameters by modeling the C IV, N V, 
and Si IV lines. During our analysis we found that the line pro- 
files in the IUE range are strongly influenced by a combination of 
three parameters at the same time: the terminal wind velocity u^,, 
the mass-loss rate M, and the impact of superionization via X-ray 
emission. It is not possible to disentangle these effects and therefore 
our solution may not be unique. 

We include the X-ray flux and temperature as deduced from 
the observations. The need to account for X-rays in wind modeling 
can be nicely demonstrated using C IV and N V doublets (Figs.|6] 
and[7](. Although it is possible to reproduce the observed C IV dou- 
blet without X-ray superionization by assuming lower M and lower 
Voa, this method fails for the N V doublet. As the effective temper- 
ature of B stars is not high enough to create sufficient amounts of 
Nv by photoionization, the observed P Cygni line profile of the 
Nv doublet can only be reproduced if superionization due to X- 
rays is included. On the other hand, the Si IV doublet in r Sco has 
little sensitivity to the X-rays (Fig.[8}. 

With fixed parameters for the X-ray emitting gas from this 
study, the terminal velocity of the wind can be deduced from fit- 
ting the line profiles. Fig.|9] shows the asymmetric line profile of 
the ClV-line doublet at AA 1548.2, 1550.8 A. We calculated three 
models taking Voo = 500, kms -1 , 1000 km s -1 and 1500km s - 
and keeping other parameters the same. Although a higher wind 
velocity cannot be excluded, we find that a model with = 
1000 km s -1 described the observed lines best. 

Setting Woo = 1000 kms -1 , attention was next directed to- 
ward a determination of the mass-loss rate. Models with mass- 
loss rates in the range log(M/MQyr -1 ) = -9.3 . . . -8.6 were 
computed and compared to the UV data. We were not able to 
find a unique solution which could describe all lines equally 
well. The C IV doublet is better reproduced with lower M = 
10 -9 ' 3 Mq yr -1 as illustrated in Fig.[l0] This mass-loss rate is also 
favored by the model fitting of the N V doublet (see Fig.ll It. On the 
other hand, Si IV doublet is best described by a model with a higher 
mass loss of M = 10" 8 - 6 Mq yr -1 as shown in Fig.H2l 

The steps as described above were performed for all our pro- 
gram stars. 



Figure 11. r Sco: Same as in Fig.[lO]but for the N V doublet. The P Cygni 
line profile of the N V doublet is strongly affected by X-ray ionization. The 
observed line is better fitted by a model with a low mass-loss rate log M = 
-9.3 (red dotted line). 
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Figure 12. r Sco: Same as in Fig.[l0]and Fig.[IT]but for the Si IV doublet. 
In contrast to the C IV and N V doublets, the observed doublet of Si IV is 
better matched by a model with a higher mass loss of log M = —8.6. 

6.3 /3Cep 

j3 Cep is a well-studied object. Despite this, to our knowledge, there 
has been no prior detai l ed mo deling of the UV resonance lines from 
this star. Do nati et al. I ( EooTh adopted a terminal wind velocity of 
l>oo = 90 km s -1 . Us ing a constraint on the product Alq(C IV) ob- 
tained bv lPrinjal j 1989b . and adopting an ionization fraction of C IV 
at 0.1%, Donati et al. estimated a mass-loss rate of log M = —9 
for j3 Cep. They also noticed that this mass-loss rate is significantly 
lower than the predicted v alue of 2.4 x 10~ 8 Mq yr -1 for a CAK 
model of a star like (3 Cep dAbbottll 1982b . 

The ionization fractions of carbon obtained by detailed mod- 
eling of the P Cep atmosphere with the PoWR code are shown in 
Fig-El The effect of ionization by X-rays prevents C IV from be- 
ing the dominant ion anywhere in the stellar wind except right at the 
photospheric level. Therefore systematically lower mass-loss rates 
would be obtained from modeling of the C IV doublet by models 
that do not account for the X-rays. (The effect is illustrated in Fig.|6] 
usin g r Sco as an exampl e.) 

ISchnerr et al. I d2008b studied wind-line variability in magnetic 
B-stars. They considered all 81 UV spectra of /3Cep available in 
the IUE archive and showed that the C IV doublet is strongly modu- 
lated with the rotation period of 12 d. Their Fig. 2 shows the gradual 
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Figure 13. Relative ionization fractions of carbon in the wind of /3 Cep 
as function of density in the wind, with declining density corresponding 
to larger radius. The wind model is calculated with Voc = 700kms — 1 , 
log M = —9.4; stellar parameters are from Table|6] 
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Figure 14. /3 Cep: Detail of the UV spectrum observed with IUE (blue thin 
line) vs. PoWR models with v^, = 700 km s -1 and log M = —9.4 (red 
dotted line). The C IV AA1548.2, 1550.8 A doublet is shown. 



transition from an enhanced to a reduced contribution of emission 
centered close to zero velocity in C IV, as typical in magnetic B 
stars. 

Figure[T4l shows the observed C IV doublet compared to a 
model that assumes Vao — 700 km s -1 and log Al = —9.4. A 
higher mass-loss rate of log M = —9.1 is required to reproduce 
Si IV line, as shown in Fig.l 151 



6.4 CMa 

This is the first analysis of the UV spectral lines of the /3 Cep-type 
star £ CMa by means of stellar atmosphere models. 

In modeling the spectrum of this object, we note that the IUE 
data of this star are of rather low quality. It seems that the back- 
ground subtraction was not performed correctly, because the mini- 
mum of the interstellar Lyman alpha absorption line is not "black" 
as it should be. Therefore we adjusted the background level ac- 
cordingly. Our model reproduces the broad SED from the UV to 
the IR very well (see Fig.|5j. Figure|4] shows our model compared 
to the optical spectrum of £ CMa - the model produces quite a 
good match to the line spectra. Therefore, it is surprising that we 
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Figure 15. ft Cep: The same as in Figs. l 141 and but for the Si IV 
AA1393.8, 1402.8 A doublet. A PoWR model with logM = -9.1 (red 
dotted line) is shown. 



2.0 



1.5 - 



■3 1.0 



0.5 




0.0 



Figure 16. g 1 CMa: Detail of the UV spectrum observed with IUE (blue 
thin line) vs. PoWR models with Voo = 700 km s -1 and logM = -10 
(red dotted line). The C IV AA1393.8, 1402.8 A doublet is shown. 



were not able to achieve good quality fits to the C IV (see Figs.l 1 6I> 
and N V lines. 

It is tempting to suggest that the reason for these model- 
ing difficulties may be due to the pole-on orientation of £ CMa 
dHubrig et al.ll201 it) , although the poor quality of the data is a con- 
cern as well. The observed line profiles of Si IV (see Fig.l 17t are 
roughly reproduced by our model. We adopted a larger turbulence 
in the wind compared to the photosphere and convolved the UV 
model spectra with a Gaussian of 1 A FWHM, corresponding to a 
turbulence velocity of about 200 km s -1 . 

For /3 Cep and V2052 Oph, the mass-loss rate required to re- 
produce the Si IV doublet is higher as compared to what is needed 
to fit C IV, an effect that we attribute to the superionization by X- 
rays. This, however, is not the case for ^ CMa. Fi gure[T6] shows 
the C IV line compared to model line with log Al = —10, as found 
fro m the analysis of Si IV . 

ISchnerr et al. I d2008l) presented a time-series analysis of IUE 
observations of the C IV doublet in £ CMa and noticed the lack of 
temporal modulations in the spectra. Increasing the mass-loss rate 
would result in a stronger absorption feature, which is not observed. 
Therefore, we conclude that the UV lines in £ x CMa are peculiar 
compared to other magnetic stars in our sample, either due to the 
quality of the observations or reflecting intrinsic peculiarity to the 
source itself. 
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Figure 17. ^CMa: The same as in Figs. 1161 but for the Si IV 
AA1393.8, 1402.8 A doublet. 
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Figure 18. V2052 Oph: Detail of the UV spectrum observed with IUE (blue 
thin line) vs. PoWR models with = 700 kms -1 and log M = -10.7 
(red dotted line). The C IV AA1548.2, 1550.8 A doublet is shown. 



6.5 V2052Oph 

This is the first analysis of wind properties of V2052 Oph. As 
with the previous stars (except r Sco), we adopt a wind termi- 
nal velocity of Voo — 700 km s -1 . We did calculate models with 
Woo = 1000 km s -1 , but results were similar. Fig. 1181 shows the 
C IV doublet in the IUE spectrum of V2052 Oph compared to a 
wind model with logM = —10.4. The effective temperature of 
V2052 Oph is relatively low at T cff = 23 kK, and C III would have 
been the leading ionization stage if there were no X-ray emission. 
As with other stars in our sample, the effect of superionization 
proves critical, in this case elevating the ionization fraction of C IV 
in the wind. 

We were not able to achieve a satisfactory fit to the N V dou- 
blet with a model that includes X-ray emission. This line is very 
sensitive to the presence of ionizing photons in the wind, and even 
models with a mass loss as low as log M = —10.7 produces a N V 
doublet that is stronger than observed. This is an interesting chal- 
lenge because the observed level of X-ray emission in V2052 Oph 
is quite low, smaller than in the other stars of our sample. 

Similarly to results for other program stars, the Si IV doublet is 
better reproduced with models that assume higher mass-loss rates 
then required for the C IV and N V doublets. Figure[T9] shows a 
model fit with log M = -9.7 to the Si IV lines. Note that the low 
T e ff of about 23 kK would normally suggest that the dominant ion 
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Figure 19. V2052Oph: The same as in Fig.QI] but for the Si IV 
AA1393.8, 1402.8 A doublet. A PoWR model with Vaa = 700 km s" 1 and 
log M = —9.7 (red dotted line) is shown. 



2.0 



1.5 - 



•3 1.0 



0.5 




0.0 



Figure 20. f Cas: Detail of the UV spectrum observed with IUE (blue thin 
line) vs. PoWR models with i>oo = 700 km s -1 and log M = —11 (red 
dotted line). The Civ AA1548.2, 1550.8 A doublet shows a purely photo- 
spheric absorption profile without any wind signature. 



stage of silicon would be Si IV; however, with the presence of X- 
rays, Si V becomes the dominant ion, and the Si IV doublet shows a 
more photospheric absorption profile. 



6.6 C Cas 

Ours is the first analysis of the wind properties of £ Cas. This star 
is the coolest among our sample. Again, we fixed the parameters 
of X-ray emission based on XMM-Newton data, and calculated a 
range of models for various values of Voo and M. 

Figure [20] shows the observed C IV doublet plotted together 
with a PoWR model. The observed line shows a purely photo- 
spheric absorption profile without any wind signature and is quite 
well fitted with a PoWR model. From such a photospheric absorp- 
tion profile, it is not possible to infer the wind velocity. We have 
therefore assumed a terminal wind speed of Voo — 700 km s _1 . 

Wind emission signatures are present in the Nv doublet as 
seen in Fig.|2T] Note that the model N V line has a P Cygni line 
profile because N V becomes the leading ionization stage via the 
effect of superionization. However, in order to reproduce this line 
we have to adopt a small mass-loss rate of only log M — —11. Sur- 
prisingly, even with this mass-loss rate, the model line is somewhat 
stronger than observed. 
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Figure 21. C Cas: Same as in Fig.|20]but for the N V AA1238.81, 1242.! 
doublet. 
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Figure 22. £ Cas: The same as in Figs.|20] and [21] but for the Si IV 
AA1393.8, 1402.8 A doublet. A PoWR model with Voo = 700 km s" 1 and 
log M = —9.7 (red dotted line) is shown. 

Similarly to other stars, there is no unique solution for the 
mass-loss rate, with different UV lines indicating somewhat dif- 
ferent mass-loss rates. Figurel22l shows models of Si IV compared 
with the IUE data. At a low T» of about 21 kK, Si IV would be the 
dominant ion stage, but X-ray emissions make Si V the dominant 
ion in the wind so that the Si IV line shows a more photospheric 
absorption profile. 



6.7 Uncertainties on M determination in non-supergiant B 
stars 

The results of the analysis of UV and optical spectra of magnetic B 
stars are summarized in the Table|6] There is a discrepancy between 
modeling the C IV and N V doublets versus the Si IV doublet, there- 
fore we determine the mass-loss rates in our program stars only to 
within a factor of a few. 

The C IV, N V, and Si IV doublets differ in their sensitivity to 
stellar atmosphere parameters. Our models show that C IV can be 
the leading ionization stage in the hottest stars in our sample, such 
as tSco and £ CMa. Then X-rays destroy C IV by photoioniza- 
tio n (see Fig.[6t. (In the context of dwarf O-sta rs this was discussed 
bv lMartins et alj|2005l : iMarcolino et al.ll2009p . In cooler stars with 
T ft < 27 kK, the leading ionization stage of carbon in the wind 
becomes C III. In this case Auger ionization by X-rays leads to the 
creation of C V, followed by recombination to C IV. Therefore, the 



C IV line is quite sensitive to both the stellar radiation field as de- 
termined by T e fi and to the X-ray emission. 

N V is not produced in the winds of stars with T» < 50 kK 
(see Fig. 3 in lHamann & Gra fener 2004]). All stars in our sample are 
significantly cooler, therefore the N V line in their winds originates 
exclusively from the ionization by X-rays. For stellar temperatures 
below 27 kK, the leading ion in the wind is N III, which is Auger 
ionized by X-rays to N V. At 27 kK and above, the leading ion is 
N IV. When X-rays are present in the wind, the leading ion becomes 
N V, or for r Sco (T» = 31 kK) even N VI. 

All models that include X-ray emission and have appreciable 
mass-loss rates show N V with a strong P Cygni line profile, which 
is not observed (see Figs.ll llandl21t.To achieve agreement with the 
observations low mass-loss rates must be adopted. For the j3 Cep- 
type variables among out program stars, the models reproduce the 
observed N V doublet best if the X-ray emission is "switched off". 
Since X-ray emission is observed, one potential explanation for cu- 
riosity could be in assuming some form of shielding of the cool 
stellar wind from the X-ray radiation. 

The ionization of Si is also sensitive to stellar and X-ray ra- 
diative field. Models predict that in the stars with T* < 27 kK, the 
leading ion would be Si IV, but the presence of X-rays make Si V 
the dominant ion. Interestingly in hotter stars with photospheric 
temperatures above 27 kK, the dominant ion is Si V with or with- 
out X-rays; the X-rays are not sufficient to produce a significant 
amount of Si VI. Let's consider ft Cep and r Sco as examples. Ig- 
noring X-ray emission in the models, Si IV is the dominant ion in 
the outer wind zone, and therefore the Si IV doublet becomes asym- 
metric in shape. With X-rays most of the Si IV in the outer wind is 
destroyed (ionized to Si V), and the remaining Si IV absoiption is 
only photospheric. On the other hand in rSco that is hotter, the 
dominant ionization stage is Si V and the Si IV line is not sensitive 
to the X-ray emission. 

These considerations show that the stellar radiation field and 
the X-ray emission have to be well known to use the UV lines as 
mass-loss rate diagnostics. There are several factors that can influ- 
ence mass-loss rate determinations, and it is worth noting these and 
their relative importance as briefly discussed below. 

Although the stellar radiation field is well described by 
our models (see Fig . ff}, the detailed p hotospheric models (e.g. 
iMorel & Butieil l2008l ; lMorel et al.l 120081) 1 give temperatures with 
uncertainties of several thousand Kelvin. Our models show that 
such differences in temperature are sufficient to alter the ionization 
stratification significantly. Thus, uncertainties in T e g are a source 
of uncertainty in M. 

For our mod els we adopted abundances as derived by 
IMorel et"ai] d2008l) . The uncertainty in abundances would reflect 
on the mass-loss rate determinations. However, we determine mass- 
loss rates only within a factor of few, which is a higher uncertainty 
than could be explained by the abundances. 

Another source of uncertainty is represented by the spectral 
and spatial specifics of the X-ray radiation. From observations we 
constrain the X-ray flux quite well; however, the X-ray spectral dis- 
tribution is only roughly known for our sources, especially for the 
X-ray fainter stars (see Table|3)- In addition, the location of the X- 
ray plasma throughout the wind is unclear. 

Similar to other non-LTE stellar atmosphere models, the 
PoWR model assumes that the statistical equilibrium is established 
locally. In a situation when the densities and the wind flow times 
are small, the timescale for recombinations can become longer than 
the timesc ale for transpor t by ad vection. In their studies of low den- 
sity winds. iMartins et ail J2005h tested the effects of advection and 
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Table 6. Stellar and wind parameters of r Sco, /3 Cep, V2052 Oph, CMa, and £ Cas. The highest and the lowest mass-loss rates M as obtained from the 
modeling of C IV and Si IV doublets for each star are given. The M value for CMa is as obtained from the model of the Si IV doublet only (see Section 5.4 
for details). The work ratio Q (see Eq.[2]for definition) is computed for the highest mass-loss rate that is given in this table. 
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adiabatic cooling on the ionization structure and the UV lines for 
a star with log Al ~ —9. They found that the inclusion of these 
effects would result in an increase of the empirically derived log 
M by ~0.15. This is within the error margin of the M derived for 
our sample stars. It is also important to notice that our final results 
are based on unclumped wind models. If the winds are clumped, 
the recombination time scale in dense clumps becomes shorter. As 
a test we calculated a grid of models that account for microlump- 
ing for j3 Cep. However, no significant difference was found com- 
pared to the un-clumped models. Moreover, in our models, the N V 
and C V produced via Auger ionization are the leading ionization 
stages, therefore our assumption of local ionization balance does 
not introduce a large error in the model ionization structure. 

We assume a spherically symmetric winds. Except r Sco that 
has a complex magnetic topology, the stars in our sample are 
oblique magnetic rotators. Their winds are likely confined by mag- 
netic fields, therefore, the treatment by the spherically symmetric 
models can be a source of uncertainty. 

Four of our pr ogram stars are j3 Cep type stars. As noticed by 
iDonati et alj ( l200ll) . in principle, the mass-loss rate may change 
by ~ 20% in response to the ~ 9% pulsation-induced luminosity 
changes, while the terminal velocity should remain roughly con- 
stant. However, these changes in the M are within the errors of our 
estimates. 

Ultimately it appears that better determinations of the effective 
photospheric temperature along with a better understanding of the 
flow structure and X-ray sources in the systems will be needed to 
obtain consistent mass-loss rates that result in self-consistent fits 
for all of the UV wind lines. 

6.8 The weak wind problem 

It is known that O-type dwarfs with luminosities log L/Z/0<5.2 
have mass-loss rates that are or ders of magnitude lowe r than pre- 
dicted by the CAK theory (see iMarcolino et al .1 d2009l ) and refer- 
ences therein). This discrepancy is commonly refereed to as "the 
wea k-wind proble m". 

iBabell ll996l) refined the CAK theory for the B-type stars by 
studying the theoretical effects of shadowing by photospheric lines 
on radiative acceleration. It was shown that these effects have large 
consequences on the winds of main sequence B stars. In particu- 
lar, the main difference from the predictions of the CAK theory 
is found for stars with T e g « 20000 - 23000 K and log g w 
3.7 — 4.0. For these objects, the mass-loss rate is found to be lower 
than predicted by CAK theory by at least a factor of 4. Among the 
stars in our sample, £Cas has values of T c g a nd log g for w hich 
the theoretical wind models were computed by |Babell ( fT996h . The 
predicted mass-loss rate for this star is 10 -8 ' 8 Mq yr _1 with a ter- 



minal wind velocity of Voo = 3.4« osc (compared to log M — 
—8.2, Voo = 1.8u sc as predicted by the CAK theory). 

Our results show that the mass-loss rate in C, Cas does not ex- 
ce ed 10~ 9 7 Mp , yr~ 1 . This is ~ 8 times smaller that the prediction 
of lBabell i fl996l) . Furthermore, we do not find any evidence for the 
predicted fast wind velocity of Voo — 2100 km s -1 . Similarly, for 
all other stars in our sample, the ma ss-loss rates a re an order of 
magnitude lower than predictions by lAbbota ( 119821) b ased on the 
CAK theory. The discrepancy with the predictions of Ivink et al.1 
d2000h is even larger. Thus we conclude that all stars in our sample 
belong to the category of weak- wind stars. Q 

This conclusion is corroborated by the high work ratio Q we 
find in our program stars (see Table|6}. As the work ratio Q we de- 
fine the mechanical work per unit time done by the radiation field 
compared to the mechanical luminosity of the wind., The exact def- 
inition of Q as computed in the PoWR code is 

= I [ffrad(r)- ^yg] dr 
Q = f r[v £ + c^ ]dr ' (2) 

where (j ra d is the radiative acceleration, P g is the gas pressure, 
and other symbols have their usual meanings. In a hydrodynami- 
cally consistent model Q = 1. When Q > 1 the model predicts 
that there is sufficient line opacity to produce radiative accelera- 
tion that is capable to drive stronger wind than the observed from 
the UV resonant doublets. We calculated hydrodynamically con- 
sistent models by choosing M values that that yield work ratios of 
Qml. Such models yield log M(Q « 1) = -9.4 for V2052 Oph 
and log M(Q »1) = -8.2 for f 1 CMa. 

Importantly, these work ratios Q are calculated f rom models 
that in clude X-ray emission in the ionization balance. iDrew et all 

1 19941) highlighted the effect that X-ray emission may have on the 
wind velocity and mass-loss rate. They suggest that ionization by 
X-rays may change the ionization structure in the inner part of 
the wind, thus reducing the total radiative acceleration and con- 
sequently the mass-loss rate. Our models include the ionization by 
X-rays in the wind acceleration zone at the level and temperatures 
indicated by the X-ray observations. Yet, we find that there is still 
sufficient radiative acceleration to drive mass-loss in excess of the 
value inferred from the UV line-fitting analysis. Therefore, the ion- 
ization by X-rays cannot be the unique solution of the weak-wind 
prob lem. 

ICassinelli et al J ( 1 1 9941) . ICassinellil d 1 994h , ICohen et all ( 1 1 997b 

2 Recently |Luc\lfeO10l) identified a weak- wind domain on a log g-log T c ff 
diagram where a star's rate of mass loss by a radiatively-driven wind is 
less than that due to nuclear burning. All our program stars belong to this 
domain. 
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studied the X-ray emission from near-main sequence B-type stars 
based on Rosat data. They found that these stars show departure 
from the ix ~ 10~ 7 Z/boi law which holds for O-type stars. From 
evaluating the emission measure of X-ray emitting gas, they con- 
cluded that a major fraction of the wind emission measure is hot, 
whereas in shocked wind theory less than 10 percent of the wind 
emission measure should be hot. 

Our observations and models support these conclusions. The 
models show that the Emission Measure 'filling factor' of the hot 
material (as determined from X-ray observations) relative to the 
cool wind Xm exceed unity for all our stars. While Xmi ~ 8 for 
r Sco is smallest, it is as large as Xgu w 200 for £ Cas and ex- 
tremely large Xmi ~ 1000 for f 1 CMa. Thus, not only are the UV 
lines in £ CMa the most difficult to reproduce from the models, its 
X-ray filling factor and wind work ratio Q are also outstandingly 
high. 

The high X-ray filling factors imply that either the hot gas 
EM^oi occupies a much larger volume than the cool wind EM coo \ 
where the UV lines are formed, or that the X-ray emitting gas has 
higher density. Both possibilities seem plausible in the case of mag- 
netic stars. The cool wind can occupy a relatively smaller volume 
because it emerges primarily from the magnetic polar regions in 
stars with dipole fields. On the other hand, high density hot plasma 
can occur in magnetically confined wind zones or loops. 



X-ray flux. The X-ray spectral temperatures of magnetic B-stars in 
our sample are not especially high. £ Cas has one of the softest X- 
ray spectra among those measured in OB stars. Moreover, X-ray 
emissions from both r Sco and /3 Cep were monitored throughout a 
rotation period, an d neither display rotational modulatio ns in their 
X-ray light curves dFavata et al . 112009k llgnace et al.ll2010h . 

Compared to non-magnetic stars, the harder X-ray spectrum 
and rotationally modulated X-ray variability are predicted by dif- 
ferent mod els of hot plasma production in m agnetically confined 
winds (e.g. Ili et alj|2008l iGagne et al1l2005l) . Can these models, 
which assume that wind motion is governed by the magnetic field, 
be applied to our program stars? 

It is usual to describe the relative importance of magnetic 
fields in gases by the plasma-/3 parameter with /3 P = Srvp/B 2 , 
where p is the gas pressure. The gas is magnetically dominated 
when P P < 1. 

For supersonic flows such as stellar winds the ram pressure 
exceeds the gas pressure and the dynamical importance of a mag- 
netic field is defined by the ratio of wind kinetic to magnetic energy 
density, given by the supersonic flow /3 viz. Airpv 2 /B 2 , where v is 
the supersonic flow speed. A small or large value indicates whether 
the magnetic field locally dominates the bulk motion, or vice versa. 
Supersonic rotation can be treated similarly (e.g. Brown et al 2008). 
Specifically lAltschuler & Newkirkl d 19691 and ref. therein) show 
that beyond the Alfven radius _Ra, where 



7 DISCUSSION 

7.1 X-ray emission from magnetic early-type B stars 

One of the results of our study is that X-ray properties of newly 
discovered magnetic B-type stars are diverse. Except for r Sco and 
CMa, our program stars are not especially X-ray luminous. On 
the contrary, V2052Oph and £Cas are intrinsically X-ray faint. 
Can these differences in luminosity be exp l ained i n the framework 
of the MCWS model? lBabel & Montmerlel dl997bh obtained a sim- 
ple scaling relation between X-ray luminosity and the magnetic and 
stellar wind parameters in the framework of their MCWS model. 
They predicted that the X-ray luminosity scales approximately with 
the product of wind momentum and some power of the magnetic 
field strength, L x tx B 0A Mv oo (see their Eqs. (10)-(1 1)). We ap- 
ply this scaling relation to predict X-ray luminosities for our sample 
stars with dipole magnetic fields. Among our sample stars, CMa 
has the strongest magnetic field, while its win d momentum is com- 
parable to other stars. The scaling relation of iBabel & Montmerlel 
(1997b) predicts rather weak dependence on the magnetic field 
strength. The predicted value of Lx (4 x 10 30 erg s _1 ) is an order of 
magnitude lower that the observed. On the other hand, V2052 Oph 
and £ Cas have a low wind kinetic energy and a weak magnetic 
field. Using the lower values of M from Table[6] the predicted X- 
ray luminosities for these stars (« few x 10 29 erg s _1 ) agrees well 
with the observed ones. Similarly, for the lower value of M for 
P Cep, the predicted and observed l uminosities agree qu i te well . 
Thus, the basic scalings obtained bv lBabel & Montmerlel (1997b) 
for the X-ray luminosity can qualitatively explain the difference in 
the level of X-ray luminosity among /3 Cep, V2052 Oph, and ( Cas. 
No te, however, that usin g the up per values for A/ for these stars in 
the IBabel & Montmerlel Jl997bl) scaling relation, would result in 
the X-ray luminosities that are too high compared to the observed 
ones. 

The more serious challenge is to understand the comparatively 
low temperatures of X-ray emitting plasma obtained from the anal- 
yses of the observed spectra and the lack of time modulations of the 



B 2 /8tt = pv 2 /2 



(3) 



is met, the radial stellar wind forces magnetic field lines to to be- 
come approximately radial. By contrast for R < Ra, the wind flow 
is confined by the magnetic field. 

In the context of stellar winds, wind confinement was consid- 
ered by Babel & Montmerle (1997), ud-Doula & Owocki (2002), 
Brown et al (2008) and others. Recalling that the wind density 
p — M I 'A-kv(R)R 2 , one can express Ra from the condition set 
by Eq. (3): 



Ra = 



Mv(R) 



(4) 



In case of a dipole magnetic field B — Bo(R t /R) 3 Eq. (O can be 
re-written as 



ta ■ w(r) 4 = 



( RlB 2 



(5) 



where r — R/R*, and the standard parametric wind velocity law 



1^ 



-, w(r) is assumed. It is convenient 



with v(r) = Vac (l—r~ ) 
to express 77, in terms of normalized stellar wind parameters, with 
M_9 the mass loss rate in units of 1O _9 M0 yr _1 , Vao in kms" 1 , 
and 1Z, the radius in units of solar radius Rq . Then 
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For stars with dipole fields (i.e. excluding r Sco) the field 
strength at the magnetic pole is shown in Table|6] To estimate 
rA we use as a characteristic value the field strength at the mag- 
netic equator (i.e. half of the polar field strength). Then we find 
77„(/3Cep)>3x 10 3 ,r;*(^ 1 CMa) fa 6 x 10 6 , ??«(V2052Oph)>4 x 
10 3 , and 77* (£ Cas)>6 x 10 3 . Consequently, the Alfven radius is at 
w 50 -R* for £ x CMa and <10 R r for other stars. These estimates 
are based on the highest mass-loss rates obtained from fitting of the 
UV lines. 
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However, we know that a significant fraction of the wind 
mass loss can be in the form of hot X-ray emitting plasma. We 
may hypothesize that the "true" total amount of matter is much 
higher. We calculated a hydrodynamically consistent model for 
V2052 Oph and £; 1 CMa (see Section RT8l > to obtain a mass-loss rate 
logM(Q » 1) = -9.4 for the former and log M(Q w 1) = -8.2 
for the latter. Even with these higher mass-loss rates, the magnetic 
fields regulate the wind motion up to ~ 7 i?» for V2052 Oph and 
« 20 R* for £ J CMa. 

Thus the winds of our sample stars are strongly confined out 
to several stellar radii where the wind velocity should reach its ter- 
minal value. Due to the very weak signatures of the winds seen 
in the UV line profiles, we can establish v^, only with a limited 
degree of accuracy. Using Vrx, « 700 km s -1 is reasonably con- 
sistent with the UV lines. Colliding wind streams at such a speed 
from two opposite directions of the magnetosphere should produce 
strong shocks, heating plasma up to ~ 20 MK. This is far above 
the maximum temperatures we infer from the spectral analysis (see 
Tabled). 

From the UV line modeling, we cannot exclude somewhat 
lower Voa for the winds of our stars, which would alleviate the ex- 
pectation of very hot temperatures, scaling as v^. However, there 
is, perhaps, a more severe problem of the differential emission 
measure distribution (DEM). The emission measure of the hottest 
plasma is « 35% of the total EM in £ x CMa and only «10% in 
j8 Cep. On the other hand, the plasma with the largest emission 
measure has a temperature of only about 1.3 MK in f CMa (see 
Tabled . This dominant contribution to the DEM from the cooler 
hot gas component is in apparent contradiction with the predictions 
of the MCWS models, where the D EM is expected to peak at the 
hottest temperature (e.g. see Fig. 8 in lNaze et al 

What could be the sour ce of wind heating i n magnetic j3 Cep- 
type stars? As pointed out bv lFavata et ai]||2009h . the relatively low 
temperature of the X-ray emitting plasma and the absence of flares 
makes heating due to magnetic reconnection in analogy to active 
cool stars rather implausible. Heating by wind shocks owing to the 
intrinsic LDI mechanism in radiatively driven winds is a possibility. 
However, the LDI needs further theoretical study in relation to its 
operation in complex geometries and the dynamics of magnetically 
confined winds. 

Another possibility could be heating b y the deposition of me- 
chanical energy related to stellar pulsations. ICassinelli et al. 
found from EUVE observations that the /3 Cep-type star /3 CMa has 
an EUV exceeding atmosphere model predictions by an order of 
magnitude. 0. They suggest that this may be related to the presence 
of heated regions near the stellar surface owing to pulsations and 
magnetic fields. If the upper atmosphere is heated mechanically by 
pulsational effects, it is likely that the X-ray source regions are also 
heated by a related mechanism. In this case, the dynamic time scale 
of the wind is smaller than the time required to establish ionization 
equilibrium. 

To investigate this possibility, we show in Fig.[23]the X-ray lu- 
minosities versus pulsation periods and amplitudes for few promi- 
nent ft Cep-type variables. Given the small sample it is difficult to 
draw any conclusions, yet there is no apparent dependence of the 
X-ray luminosity on period or pulsational amplitude. 



3 It is possible that this problem can be resolved by a better description of 
X-ray emission in model atmospheres (Hamann et al., in prep.). 
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Figure 23. X-ray luminosity of some prominent ft Cep-type stars vs. mag- 
nitude of puls ation (upper panel) and pul sational period (lower panel) from 
the catalog of lStankov & Handler! feOOSl) . 



7.2 Comparison of X-ray properties between B stars with 
and without confirmed magnetic fields 

In this section we compare X-ray properties of "normal" early-type 
B stars where magnetic fields have not been detected up to now. The 
XMM-Newton spectrum of the /3 Cep-type variable j3 Cen (B1III) is 
well fitted with a three temperature spectral model, with the hottest 
comp o nent fcr max sa 0.6 k eV and (kT) = 0.3 keV i Raassen et al.l 
120051) . iHubrig et al.l d2006l) and lSilvester et al.l feOQsl) searched for 
a magnetic field in the (3 Cep-type /3CMa (B1II-III) and did not 
detect magnetism, with errors on a longitudal field of just 10 G. 
j3 CMa was observed by XMM-Newton, its spectrum is well de- 
scribed by a multi-temperature model with the hottest tempera- 
ture component at kT = 0.3 keV (Waldron et al. 2011, in prep.). 
Chandra LETGS spec trum of aVir (B1III-IV) was analyzed in 
IZhekov & Pallal d2007h . who report the maximum of the emission 
measure distribution as being between 0.2 keV and 0.3 keV. We ob- 
tained LETGS Chandra observations of aCru (B0.5IV). Its X-ray 
temperature and luminosity are similar to those of aVir (Ignace 
et al. in prep.). j 3 Cm is another Bl giant observed by Chandra. 
IZhekov & Pallal J2007h find the maximum of the emission mea- 
sure distribution in this /3 Cep-type variable is between 0.2 keV and 
0.3 keV. Thus, it appears that X-ray luminosities and temperatures 
are very similar among the /3 Cep-type stars. 

On the other hand, the X-ray luminosities of magnetic stars 
have much greater dispersion, with two significant outliers: the X- 
ray luminous £ CMa and the X-ray faint V2052Oph, while their 
spectral X-ray temperatures are rather similar (see Table|3}. Due to 
the small number of stars included in our investigation, this conclu- 
sion is only preliminary. Larger samples need to be considered to 
draw firm conclusions about the differences in X-ray properties of 
magnetic and "normal" early-type B stars. 

Overall, it has emerged from our analysis that high X-ray lu- 
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minosity and hard X-ray spectra are not necessary observational 
characteristics of the magnetic early-type stars. This conclusion 
fully confirms the earlier results based on t he RASS data (see 
also r ecent results on X-rays from AOp stars dRobrade & Scnmittl 
l201lk 

One of the complications in massive early-type star studies 
is that a significant fraction of such stars is in binary or multiple 
systems and may have spatially unresolved companions that are in- 
trinsic X-ray sources. An X-ray source spatially coincident with a 
B star can be due to a lower m ass companion dCzesla & SchmitJ 
l2007l : lRemage Evans et alfcOllh . In this paper we consider B0-B2 
type stars. These stars are young, not older that a few 10 Myr. Possi- 
ble low-mass companions would be c oronally active with X-r ay lu- 
minosities of 10 28 -10 30 ergs" 1 (e.g. lFlaccomio et alj|2003f) . Two 
newly X-ray detected stars considered in this paper, V2052 Oph and 
( Cas, have X-ray luminosities in this range (see Table[4). There- 
fore, in principle, we cannot rule a confusion with an unseen coro- 
nal component. Nevertheless, we have strong indications that the 
observed X-ray emission is intrinsic for the early-type B stars. The 
spectra of X-ray emission from £ Cas and V2052 Oph are soft, e.g. 
in £ Cas the average temperature is 1 MK. This is not usual for 
young TTau typ e stars, which have average temperatures in the 
range 5-30 MK dGiidel & Nazj|2009l) . The temperature of the X- 
ray emitting plasma in V2052 Oph is similar to other f3 Cep-type 
stars, albeit its X-ray luminosity is quite low (see Fig.123). From 
the analysis of the UV spectra of our program stars, we know that 
X-rays are present in the winds since they are required to enhance 
the fraction of C IV and N V in the winds (see Sect. l6.5t . On this 
basis, we believe that the observed X-ray emission is intrinsic for 
/3Cep, V2052Oph, (Cas, and CMa. Among stars in our sam- 
ple, j3 Cep is a known binary star with a late Be-type companion. 
Late Be stars are not known to be intrinsic X-ray sources, therefore 
it seems reasonable to assume that the observed X-ray emission is 
intrinsic to the B2 component. 

7.3 X-rays in stars with very small M 

C, Cas displays quite a soft X-ray spectrum with nearly 85% of the 
EM at temperature <1 MK (Fig.[3]and Table[3]). Our estimate of its 
mass-loss rate (M ~ 10" 11 M0yr -1 ) places this object in the 
regi me where the winds may d i splay u nusual properties. 

ISpringmann & Pauldrachl d 19921) demonstrated that the as- 
sumption of a one-component fluid will not be valid for low den- 
sity winds. Instead, the metal ions lose their dynamical coupling to 
the ions of hydrogen and helium. The metal ions will move with 
high velocities, with helium and hydrogen failing to be dragged 
along. The collisionally induced momentum transfer is accompa- 
nied by fri ctional heating, which dominate s the energy balance. 
As a result, ISpringmann & Pauldrachl d 19921) predict electron tem- 
peratures of T ~ 1 MK in the outer wind regions for stars with 
m ass-loss rates of ~ ^q- s...9 yr -1 . The subseq uent analysis 
o flKrticka & Kubad feOOCl) and lOwocki & Pulj i2002T) favored even 
weaker winds with M ~ 10 _11 A/q yr _1 for the wind decoupling 
to occur. 

The theoretica l analysis of this p r oblem is not settled (e.g. 
Votru ba et al ] |2007h . lGavlev & Owockil i ll 9941) investigated the ef- 
fects of "Doppler heating" in stars with low mass-los s rates , 
and confirmed the conclusions of lSpringmann & PauldracrJdl992h . 
They also noticed that Doppler heating may lead to an instability, 
possibly resulting in a larger degree of microturbulence. This is an 
interesting prediction since our analysis favors large microturbulent 
broadening, as in the case of i 1 CMa. 



Given its small M value and very soft X-ray spectrum, it is 
tempting to suggest that C, Cas may be a candidate object for wind 
decoupling. 

Another interesting mechanism for gener ating X-ray emissio n 
in rotating magnetic stars was put forward bv lSuzuki et al.M2006h . 
They proposed that extended X-ray-emitting regions can exist in 
massive stars with thin winds owing to collisionless damping of 
fast MHD waves. Stellar rotation causes magnetic field lines an- 
chored at the stellar surface to form a spiral pattern, and magneto- 
rotational winds can be driven (e.g., Weber & Davis 1967). If the 
structure is magnetically dominated, fast MHD waves generated at 
the surface can propagate almost radially outward and cross the 
field lines. The propagating waves undergo collisionless damping 
owing to interactions with particles surfing on magnetic mirrors 
that are formed by the waves themselves. The dissipation of the 
wave energy produc es heating and acceleration of the outflow. The 
ISuzuki et alj (2006J) mechanism works eff ectively in moderate and 
fast rotating stars. Considering Figure 6 in lSuzuki et al.ld2006l) . the 
magnetic field strength, wind density, and rotation speed in £ Cas 
are such that the star may be in the domain where the Suzuki et al. 
mechanism is operational. 



7.4 On the multi-phase structure of winds from magnetic 
early B-type stars 

We try to understand whether some of the facts that have emerged 
from our analysis can be qualitatively explained by a picture of 
an oblique magnetic rotator where the bulk of the stellar wind is 
governed by the dipole magnetic field. 

Despite relatively strong surface fields, oblique magnetic rota- 
tors have open field regions near their magnetic poles, and thus may 
have sectors of wind flow that are similar to normal B type stars. 
The UV spectra could originate in these cool sectors of the wind 
with a smaller filling factor as compared to the hot gas. The wind 
distribution is clearly asymmetric (not even axisymmetric owing to 
rotation). Perhaps this could explain the strange properties of the 
UV lines seen in i 1 CMa that we view nearly rotational pole-on to 
the rotation axis, and possibly the large dispersion in X-ray lumi- 
nosities among our sample stars as well. 

The hot part of the wind is apparently heated up to a few MK. 
This matter is confined by the magnetic field up to 10 i?*, so it oc- 
cupies a large volume, and, in parts, may have high density. This 
would explain the large filling factors of hot material we infer from 
our models. The spatial separation between hot and cool gas may 
explain how the observed N V doublet can be weaker than predicted 
by our stellar atmosphere models. The total amount of matter in the 
cool and hot wind components may be close to values predicted by 
our hydrodynamically consistent models, and so the "true" mass- 
loss rates are a factor of a few higher than obtained from the empiri- 
cal fits of UV lines using a model that relies on spherical symmetry. 



8 SUMMARY OF RESULTS 

Dedicated XMM-Newton observations were obtained for three 
early-type magnetic B-type stars, t; 1 CMa, V2052Oph, and £Cas. 
We report first detections of X-ray emission from V2052 Oph and 
£Cas. The observations show that the low-resolution X-ray spec- 
tra of our program stars are well described by a multi-temperature 
CIE plasma. The bulk of the emission measure originates from the 
plasma with a temperature of ~ 1 MK. The X-ray luminosities dif- 
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fer by large factors, with £ CMa being the most X-ray luminous 
star in our sample, and V2052 Oph the least luminous. 

We compare X-ray properties of /3Cep-type variables that have 
magnetic field detections against those without such detections. 
Our comparison shows indications that the general X-ray proper- 
ties are quite similar, although the magnetic stars display a greater 
dispersion in their X-ray luminosities. We searched for correlations 
between X-ray emission and pulsational and rotational properties 
of the stars but were not able to find any. Larger samples need to be 
studied to draw firm conclusions. 

The basic MCWS model scaling relation between X-ray lu- 
minosity and stellar mag netic field and wind parameters 
l lBabel & Montmerlg|l997bl) qualitatively describes the difference 
in the observed level of X-ray emission from our sample stars. 

We compiled a sample of peculiar early B-type magnetic stars and 
considered their X-ray properties based on archival data. The com- 
parison shows that these stars, despite very similar stellar and mag- 
netic field parameters, have quite different X-ray properties: while 
some stars are hard and luminous X-ray emitters, others are ap- 
parently soft and rather faint. This new data confirm earlier results 
based on studies of Ap-Bp stars using RASS data. 

We analyzed spectra of five non-supergiant B stars with magnetic 
fields by means of non-LTE comprehensive stellar atmosphere code 
PoWR. The PoWR models accurately fit the stellar photospheric 
spectra and also accurately reproduce the SED from the UV to the 
IR band. 

The PoWR code was used to empirically obtain wind parameters 
of r Sco, p Cep, f 1 CMa, V2052 Oph, and ( Cas from the analysis 
of C IV, N V, and Si IV doublets in stellar UV spectra. 

The model UV lines depend sensitively on the parameters of X- 
ray emission, wind velocity, and mass-loss rate. X-ray properties 
as derived from the observed spectra were incorporated into the 
wind models. We confirm that the emission measure filling factors 
of X-ray emitting material are high, exceeding unity for all stars. 

Our analysis revealed the weak wind problem for the magnetic 
early type B stars. The wind terminal speeds are on the order of 
stellar escape speed. The inferred mass-loss rates from the UV line 
analyses are significantly lower than the theoretically expected and 
predicted by hydrodynamically consistent wind models. 

Although X-rays strongly impact the ionization structure in the 
wind, their effect does not reduce the total radiative acceleration 
enough to explain the low mass-loss rates deduced from modeling 
the UV lines. When X-rays at the observed level and temperatures 
are included in the model, there is still sufficient radiative acceler- 
ation to drive a stronger mass-loss than the observed one. 
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